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I. INTRODUCTION 

The system comprised of an electron spin coupled to a large number of nuclear spins is 
a popular object of radio spectroscopy. The relevant mathematical model related to high 
magnetic fields typical for the EPR and NMR spectroscopy is studied fairly well. The 
recent interest to this model with weak magnetic fields is associated with appearance of new 
physical objects (quantum dots) whose spin dynamics may be described, with a sufficient 
accuracy, by this model. The up-to-date laser technic makes it possible to prepare the 
electron spin in a quantum dot in the state with a specified projection, while observation of 
the polarized luminescence allows one to detect the electron spin dynamics. General picture 
of the electron spin dynamics in these experiments is described in sufficient detail in the 
appropriate theoretical papers (@, 0, H, @, 0, @, @, EH)- The spin Hamiltonian used 
in these papers is well known in spectroscopy of free radicals and similar systems. At the 
moment of our first publication of this paper in arXiv, we were not aware of any papers 
that describe and use the fact that, in the case of uniform electron density on nuclei, such a 
Hamiltonian can be diagonalized exactly. This is why, in the first version of this publication, 
we presented solution of this problem. However, just a few days after our publication, we 
received a reference to the paper [1| in which this result was described. Thus, the results 
presented below in section 2, in fact, have been already obtained earlier in Note also 
that the system treated below is characterized by a large number of the degrees of freedom, 
and exact solution of the relevant kinetic problems may be of independent interest. 

The system under consideration is comprised of electron spin 1/2 coupled to even number 
2N of nuclei with spins 1/2. Interaction of the electron is supposed to be the same for all 
the nuclei. The Hamiltonian of the problem has the form: 

H = u)S z + A\\S z I z + A x (S+I-. + S-.I+) (1) 

2N o 

I=£la, ll = J (2) 

a=l 4 

Hear S z I z -are the ^-projections of the electron spin and total nuclear angular momen- 
tum, S(I)± = S(I) X ± tS(I) y - are the standard lowering and raising operators for the z- 
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projection of the electron ( total nuclear) spin, u — is the external magnetic field expressed 
in frequency units, An A±. - are the components of the hyperfine coupling tensor. The 
main results of the paper are: 

1. It is shown that Hamiltonian (1) can be diagonalized exactly. 

2. The formulas that describe dynamics of the initially polarized electron spin for a given 
initial nuclear density matrix are derived. Note again that this result was presented in[l|]. 




3. Formulas that describe pulse-to-pulse transformation of the nuclear density matrix for 
periodic orientation of electron spin by optical pulses, which is typical for experiments with 
quantum dots, are obtained. 

4. The calculated steady-state nuclear density matrices for the regimes of constant-sign 
and sign- alternating optical orientation of the electron spin. It is shown that, in the case of 
constant-sign orientation, a nuclear polarization arises whose ultimate value is calculated. 

5. It is shown that typical experimental regimes of the electron spin optical orientation 
may be accompanied by a certain ordering of the nuclear density matrix with respect to the 
total nuclear angular momentum and its projection. The latter is reveal in quasi-periodicity 
of the electron spin dynamics, which is absent in the case of the high-temperature nuclear 
density matrix. 

6. The examples of calculations of the electron spin dynamics and magnetic- field depen- 
dence of its steady-state polarization for various regimes of optical orientation are presented. 

7. The echo effect in the dynamics of the electron spin in a quantum dot after a 7r-pulse 
of the magnetic field is predicted. 



Hamiltonian (1) depends on the total angular momentum of all the nuclei. So, it makes 
sense to employ the representation of wave functions that are eigen functions the operators 
of total nuclear angular momentum squared and its z-projection. 

These functions have the following quantum numbers: 

I - is the total angular momentum quantum number related to the mean value of the 
angular momentum squared by the formula: (I 2 ) = 1(1 + 1) (in the units of h 2 ), 

L - z-projection of the total nuclear angular momentum (in the units of h). 

This set is, however, not complete because one and the same value of the total nuclear an- 
gular momentum can be obtained by summing up the elementary nuclear angular momenta 
in many different ways. Let us number the set of the remaining quantum numbers by the 
index a: ^iL,a = \I,L,a). In this representation the matrices of the total nuclear angular 
momentum squared and its ^-projection are diagonal. Note that for the states with a given 
quantum number of the total nuclear angular momentum squared /, the diagonal elements 
of its matrix are 1(1 + 1) while the diagonal elements of I z are L = —1, 1 — I, — 1, /. 
Thus the dimesionality of the subspace with given number / and a equals 21+1. As is 
known from the theory of angular momentum, the matrices of the x- and y- projections of 
the total angular momentum I x and I y (and, therefore, I± ) have no nonzero matrix ele- 
ments between the states with non equal quantum number of the total angular momentum: 
(I, ...\I XyV \I' , ...) = for I 7^ V . It follows from the fact that the operator of the total angular 
momentum squared commutes with the operator of any of its projection. Let us show that 
the set of states \I, L, a) can be chosen in such a way that the operators I± will have zero 
matrix elements also between the states differing only by the quantum number a even for 
the same /. Consider, for this purpose the states with a given total angular momentum / 




II. DIAGONALIZATION OF THE HAMILTONIAN 
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and the maximum z-projection L = I at different a: \I, I, a). The number of these states 
will be equal to the number of ways that allow one to obtain the total angular momentum 
/ by summing up 2N angular momenta 1/2. These states are linearly independent and, 
without the loss of generality, they can be regarded as normalized and mutually orthogonal: 
(/, I, a\I, I, a') = 5 aa /.. We obtain all the rest states from these ones using subsequently the 
lowering operator for z-projection: 

\I,L,a) = r\I,I,a), L = I-n n = 0,l,...,2J (3) 

Let us show that the set of functions obtained in this way is orthogonal. Each of these 
functions, in conformity with the theory of angular momentum, is the eigen function of the 
operator of the total angular momentum z-projection with the eigen value L. For this reason, 
the functions with different quantum numbers of the total angular momentum z-projection 
will be necessarily orthogonal as eigen functions of a self-conjugate operator corresponding 
to different eigen numbers. What is left is to show orthogonality of the functions differing 
only by the quantum number a. So, we have: 



(1,1 -n,a\I,I -n,a') = (I, I, a\IlF\I, /, a') = (I, I, a\ I^I+I^V, I, a') (4) 



In conformity with the theory of angular momentum, the function J™ _1 | J, J, a') will be 
necessarily eigen function for the operators of the total angular momentum squared and its 
projection, with the eigen numbers being equal, respectively, to /(/ +1) L — I — n + 1. 

As a result, taking into account the identity 



I 2 -I +/, 



we have 



n-l 



I, I, a') 



I(I + 1)-L 2 + L 



L=I-n+l 



and, therefore, 



(/, I — n, a\I, I-n, a') ~ (/, /, a^ 1 ^ 1 ^, /, a) 
Repeating this reasoning n — 1 times more, we obtain: 



(5) 
(6) 

(7) 



(J, I — n,a\I,I — n, a') ~ (I, I, a\I, I, a') = 5 aoi i (8) 

Thus the set of functions (3) is orthogonal. Then, it follows from the above calculations 
that the operators I±, for this set of functions, have no nonzero matrix elements between the 
functions with different a. Vanishing of the matrix element of the type (/, L, a\I±\I, L=pl, a') 
for a/a' can be proved by the procedure similar to Eqs. (4) - (8). At a = a' the values 
of these matrix elements are known from the theory of angular momentum 

(I, L, a\I+\I, L - 1, a) = (I, L - 1, L, a) = yjl(l + 1) - L(L - 1) (9) 

Let us examine the matrix of Hamiltonian (1) in the representation of the functions 
obtained by multiplying the functions (3) (corresponding to the nuclear degrees of freedom) 
by the wave functions of electron spin with a given z-projection \S), S = ±1/2: \S, I, L, a) = 
\S)\I, L, a). First, as it follows from the above results, this matrix can be divided into blocks 
with different total nuclear angular momenta /. The dimension of each of the blocks is given 
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by: (the number of possible projections of the total nuclear angular momentum) x (the 
number of possible projections of the electron spin)= (21 + 1)2. The number of these blocks 
with a given total angular momentum of the nuclei is equal to the number of ways how 
the angular momentum I can be obtained by summarizing 2N elementary nuclear angular 
momenta 1/2. We denote this number by Tn(I). It can be shown (the idea how one can do 
this is presented in (3]) that 



W) = C 2 V - C 2 V _1 1 = 0,1,..., AT. (10) 

Thus to diagonalize Hamiltonian (1) it suffices to diagonalize each block with a given nuclear 
angular momentum I. Let us order the basis functions of a chosen block in the following 
way: |+, J), | — , I), |+, I — 1), ...| — , 1 — I), |+, — I), | — , — I). One can easily see that, for the 
basis functions arranged in this way, the considered block of the Hamiltonian appears to be, 
in turn, divided into blocks, with two of them, corresponding to the states |±, ±J), being 
one dimensional and 21 of them, two dimesional. General form of the two dimesional block 
corresponding to the states | — , L) |+, L — 1), where L = I, I — 1, 1 — I, is 



H LI = | + W A ±y /I(I + l)-L(L-l) f u { , i 

" A ± Jl(I + l)-L(L-l) (w + A,i(L-l))/2 



(11) 

Diagonalization of the two dimensional blocks can be easily performed and, as a result, 
for the final form of the spectrum of Hamiltonian (1) we have: 

1 / An 



E^ = -i K -^±^[u + A ll L-A ll /2f + AAl[I(I+l)-L(L-l)]j, L = I,...,1-I (12) 

, AnI±U 

E t = — 2 — > / = 0,l,...,iV 
The multiplicity of degeneracy of each eigen-energy is given by Eq. (10). 



III. DYNAMICS OF THE ELECTRON SPIN 

As was mentioned in introduction, Hamiltonian (1) can describe spin dynamics of a 
quantum dot. In the simplest way, the typical experiment can be represented as follows. 
A short polarized laser pulse transverse the quantum dot into the excited state with the 
spin-oriented excited electron. Then, the spin dynamics of this electron is detected by 
time dependence of polarization of its luminescence. This dynamics is considered to be 
determined by the contact interaction of the electron spin with the nuclei of the quantum 
dot. Hamiltonian (1) corresponds to the assumption that the excited electron wave function 
uniformly covers all the nuclei of the quantum dot. Bearing in mind this experimental 
arrangement let us assume that at the time moments t = T m ,m = 0, 1,2,3, ... the system 
is subjected to the action of short high-power optical pulses which orient the electron spin. 
Now the problem can be mathematically formulated in the following way: 

1. Right after the pulse, the electronic density matrix is assumed to have the form: 

p e (T m + 0) = Sp /P (T m + 0)= ( COS 2< ^ s . n ° 2 J 
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where the symbol Sp/ denotes diagonal summation of the total density matrix before the 
pulse p(T m + 0) over the nuclear quantum numbers. The angle <fi quantifies the degree of 
orientation of the electron spin by the optical pulse and depends on the pulse polarization. 
Thus the observables associated with the electron spin at t — T m , generally, experience a 
jump. 

2. We assume that all the observables associated with the nuclear system are, on the 
contrary, continuous at t — T m , i.e., that they are not directly affected by the optical pulse. 
It means that the nuclear density matrix 

P n = SP S P (13) 

is continuous at t = T m , i.e., p n {T m - 0) = p n {T m + 0) = p n {T m ). The symbol Sp s in (13) 
denotes diagonal summation of the total density matrix over electronic quantum numbers. 

3. The above requirements imposed on the electronic and nuclear density matrices are, 
generally, not enough to unambiguously reconstruct the total density matrix of the system 
after the pulse. To do that, one has to consider the dynamics of the system during the pulse. 
This, however, lies outside the scope of this paper. Below, we make the simplest assumption 
compatible with the above requirements that the total density matrix immediately after the 
pulse can be represented as a product of the nuclear density matrix before the pulse and the 
electronic density matrix corresponding to the oriented electron spin: 

P (T m + )= P " { T m )(™l* J, J (14) 

4. We assume that at the very beginning of the experiment the nuclear density matrix 
corresponds to equi-probable distribution over all nuclear states (high-temperature density 
matrix) . 

5. We seek for the dynamics of the system under the action of a train of optical pulses. 
The problem formulated in this way can be consistently solved. First of all we will show 

that the nuclear density matrix, under the assumptions made above, is always diagonal 
in the representation of functions with specified total nuclear angular momentum and its 
projection. Let at t — 0, the system to be excited by the first optical pulse that orients the 
electron spin. Then, at t — +0 the electron spin appears to be 1/2 (cos 2 — sin 2 0), while the 
nuclear state remains unperturbed (high-temperature state). Let us write the appropriate 
density matrix (14) in the representation of the functions \S,I,L,a) used in the previous 
section: 

(SJ,L,a\p\S , J\L\ar} = 2-™6ss>6ir6 w 6 c «{^ 2 * gI-l/2 (15) 

In accordance with the above results Hamiltonian (1) can be divided into blocks of total 
nuclear angular momenta I (referred to as /-blocks), with two one dimensional and 21 two- 
dimensional blocks within each of them. We will call last-mentioned two dimensional blocks 
L/-blocks, with L being the projection of the total nuclear angular momentum corresponding 
to the upper-left corner of the block (11). To have a clearer idea about the topography of the 
matrices considered in the problem under study, we show in fig. 1 the block corresponding 
to the states with 1 = 1. The symbols " + " and " — " denote nonzero elements of the 
density matrix (15) in this block, related to the electron spin projections +1/2 and —1/2, 
respectively. The nonzero elements of Hamiltonian (1) are colored gray. Since the initial 
density matrix (15) is diagonal, the nonzero elements of the density matrix may appear, in 
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process of its evolution only where the matrix elements of the Hamiltonian are nonzero, i.e, 
in the places colored gray. Thus the total density matrix under the initial condition (15), 
will have the topography identical to that of the Hamiltonian, i.e, it will be block- diagonal. 
Now let us assume that the system under consideration is acted upon by the second optical 
pulse polarizing the electron spin. According to the formulation of the problem described 
above the density matrix of the system will have the form Eq.(14). To obtain the element 
p\ v of the nuclear density matrix, entering Eq.(14), we have to sum up the elements of the 
total density matrix on the diagonal of the two dimensional block corresponding to these 
nuclear states. For instance to obtain the matrix element p\_\ (fig-1) we have to perform 
diagonal summation in the upper-right block (the relevant matrix elements of the total 
density matrix are located in the crossed cells) All the blocks whose diagonal summation is 
needed to obtain the nuclear density matrix are separated in fig.l by thick lines. Since the 
nonzero matrix elements of the total density matrix are located in the cells colored gray, 
one can see from the fig.l that the nuclear density matrix after the second pulse will be 
diagonal in the representation of the total nuclear angular momentum and its projection. 
Using similar reasoning to the third, forth etc. pulses, we obtain (by induction), the nuclear 
density matrix, for the problem under consideration, is always diagonal in the representation 
of the total nuclear angular momentum and its projection. Since the hole problem brakes 
down in to blocks of the total nuclear angular momentum I we can denote the diagonal 
elements of the nuclear density matrix of the block I by Uli,L = —1, 1 — I, — 1,1. 
To completely specify the block, one should have indicate its quantum number a of the 
used representation. However, we shall not present it because no observable depends on this 
number. Inside the block with a given J, one can consider the two-dimensional L/-blocks 
corresponding to those of Hamiltonian Hh (11). The corresponding Li-blocks of the density 
matrix, which are further denoted by Qli, meet the equation 

ie L i = [H L i,e LI ] (16) 

The quantum-mechanical problems related to two dimensional matrices are well known 
and, in what follows, we use the appropriate terminology. It is convenient to express the 
matrix Hu in terms of the Pauli matrices cr x ,y,z' 

H L i = — + b LI (T z + c LI a x (17) 

where 

b LI = A\ ] (l/2-L)-u c LI = 2A ±y /l(I+l)-L(L-l) (18) 

The sought density matrix Q^j can be also expressed in terms of the unit matrix (with 
the coefficient denoted by 9q 1 ) and the Pauli matrices 

Qli = 0^ + 0^cr x + 6^a y + 6?<T t (19) 

Here, the quantities 0^ y z are the components of so the so-called Bloch vector (l3| . 
Then, the standard solution of the problem of dynamics of a two-level system yields for 
the Bloch vector the following expression: 

e LI (t) = M LI (t)e LI (o) (20) 
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Where the matrix M L/ is defined as: 

(1 + gr(cosfit- 1) sin fit ^(1- cos fit) \ 



M LI (t) 



^ sin fit cos fit — ^ sin fit 

V ^(1-cosfit) ^sinfit 1 + ^(cosfit - 1)7 



(21) 



'LIi 



c LI , fi = fi L/ = V& 2 + c 2 



Let the optical pulse create the electron in a partially polarized state with the density 
matrix 



cos 2 
sin 2 



(22) 



At = (tt/2), the value of the electron spin is 1/2 (—1/2). 

For the initial state of the type (14), where the nuclear density matrix has only diagonal 
elements Uli, the initial condition for the block LI of the total density matrix has the form: 

»»(°) = ( u "f* tfJU) (*» 

9 LI : For the initial values of the Bloch vector and the quantity 9 related to the 
block LI, it yields: 

LI (0) =| ] , 9% 1 = hu LI sin 2 + U L -x,i cos 2 0) (24) 

U L i sin 2 - U L -i,i cos 2 / 

Contribution of the block LI to the mean value of the electron spin is — Sp Ou(£)cr z . 
Its convenient to use the double value of this contribution, which is denoted as 

A LI (t) = -2Sp Q LI (t)a z = -e L J (25) 



It follows from Eqs.(20),(21) and (24) that 

A LJ (t) = (V L _ 1;J cos 2 - C/ L7 sin 2 0^ (l + Jr(cos fi L/ t - 1)) (26) 

When calculating contribution of the block with a given total angular momentum /, we 
have to take into account the two one-dimensional blocks, which provide the following time- 
independent contribution to the mean value of the electron spin: (Ujj cos 2 0— sin 2 0)/2. 
Then the general formula for the time dynamics of ^-projection of the electron spin prepared 
, at t = 0, in the state (22) for the initial density matrix Uli is given by: 

1 N f 1 \ 

(S(t)) = -Y,r N (I)( £ A L i(t) + Uij cos 2 sin 2 (27) 

1 1=0 V L=1-/ ' 



A L j(t) = Ul-i,i cos 2 - Uli sin 2 



^-(cosfi L/ t- 1) 



b L i = A,|(l/2 - L) - u, cli = 2A ±v //(/ + l)-L(L-l), 



fi 



LI 



6 LJ + C LI 
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F ( J\ — siN-I—l 
1 N \ A ) — Lj 2N ~ Lj 2N 

It is often interesting to know the steady state value of the electron spin polarization. 
This value can be obtained by removing the time-dependent component in each of the 
contributions (27), i.e., using the substitution 

A LI -> = (u L - X j cos 2 <\> - U LI sin 2 (l - = (Ul-v cos 2 <\> - U LI sin 2 ^ 

(28) 

For the case of (i) the isotropic contact interaction 2A± = An = A, (ii) zero magnetic field 
to — 0, (iii) circularly polarized optical pulse = 0, (iv) high-temperature initial nuclear 
density matrix Uli = 2~ 2N , and (v) 2A ^> 1, for the electron spin dynamics, one can derive 
a simple formula. 

To do that, note that: 

1. the frequencies Qli — AI, in this case, do not depend on the quantum number of 
projection of the total angular momentum, 

2. at 2N > 1 we may assume that c 2 LI = A 2 (I 2 — L 2 ), 

3. for r N (I) given by (10), we can use the asymptotic expression 

2 2N d / I 2 \ 

Then, by substituting in Eq.(27) the summation by integration and calculating the inte- 
grals, we obtain 

1 1/ NA 2 t 2 \ ( NA 2 t 2 \ , . 

A formula similar to (30) was derived based on intuitive considerations about the nuclear 
field and its statistics in |^|. As seen from Eq.(27), the electron spin dynamics is controlled 
by the nuclear state of the system, which, in turn, can vary under the repetitive orientation 
of the electron spin by optical pulses. So in the next Section we will derive the laws of 
pulse-to-pulse transformation of the nuclear density matrix and will calculate the steady- 
state nuclear density matrices for different regimes of the repetitive optical orientation of 
the electron spin. 



IV. NUCLEAR DYNAMICS 

A. Transformation of the nuclear density matrix under optical orientation 

As was shown in the previous section, the density matrix blocks corresponding to states 
with a given total angular momentum / are mutually independent. So, in what follows we 
will consider only one of the blocks, for which we formulate the following problem: 

1. Let our system be repetitively (with the period T) acted upon by an optical pulse 
which orients the electron spin creating the states of the type (14) at the moments T m = mT. 

2. Let the block of the nuclear density matrix corresponding to the states with the total 
angular momentum I have to the moment of arrival of the n-th pulse, the diagonal elements 

jjn 
° LI- 

3. Find the elements of this block U^f to the moment of arrival of the (n + l)-th pulse. 



9 



As was already mentioned, the nuclear density matrix, in our case, is always diagonal, and 
to obtain its elements, one has to sum up diagonal elements of the diagonal two-dimensional 
block with a given value of the total nuclear angular momentum projection and with different 
projections of the electron spin. Note that these blocks do not coincide with the blocks LI. 
In fig.l, these blocks are separated by thick lines, while the blocks LI are colored gray. 
Keeping this in mind, we can express the diagonal elements of the nuclear density matrix 
in terms of the elements of L/-blocks of the density matrix and the parameters defined by 
(19): 

U u = ©22 +1,/ + ©n = 0o +M + Oo 1 ~ ~ e L J) (31) 

For the extreme values L = ±1, we have to assume that 

e£ 1>J = U H n 7 ' 7 = U- n (32) 
Using the results of the previous section and introducing the quantities 



= I ( CLI 



cos(Q LI T) - 1 



(33) 



2\n LI/ 

we obtain the following formulas for the pulse-to-pulse transformation of the nuclear 
density matrix: 

U n L j l = U n LI [l + $ LJ sin 2 + $ L+V cos 2 0] - (34) 

-Ul-i,i*li cos 2 - U2 +hI $ L+ i,i sin 2 + U ° ~ ^ L — I — l, 1 — 7 

For the the extreme values L = ±1: 

U H +1 = U H + $ 77 (U n sin 2 - j cos 2 0) + -5 ^ (35) 

r 

U — U n 

[/«+/ = [/^ _ ^^^(^ , sin 2 - ^ cos 2 0) + -2 =^ (36) 

r 

The last terms in Eqs. (34) -(36) take into account the possible nuclear relaxation, with 
r being the relaxation time in the units of T and U = 2~ 2N being the diagonal elements of 
the high-temperature nuclear density matrix. 



B. The steady-state nuclear density matrix for the case of repetitive sign-constant 

orientation of the electron spin. 

If the angle does not change from pulse to pulse, then after sufficiently long time, the 
nuclear density matrix will stabilize lim^oo U^j = Ufj. The equations for the steady-state 
nuclear density matrix can be obtain from the conditions U^j 1 = = C/£j 

U* [® LI sin 2 + $ L+M cos 2 0] - Ul\^ LI cos 2 - U? +1J <S> W sin 2 + U ° - U ™ = (37) 



L = I I 
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For the extreme values, L = ±1: 

sin 2 - Uf_ x j cos 2 0) + U ° - U " = (38) 



T 



- Qi-ijiU^u sin 2 - C/** 7 7 cos 2 0) + U ° = (39) 

Since the blocks / are independent the condition of normalization 

]T U LI = (21 + 1)U I = 0, 1, 2, N (40) 

L=-7 

should be met. Normalization of the total nuclear density matrix has the form: 

N I 

E r iv( 7 ) E u LI = i (4i) 

7=0 L=-I 

It can be easily verified that in the absence of relaxation r — > oo and = (0 = vr/2), 
the steady state nuclear density matrix corresponds to the complete polarization of the 
nuclear angular momentum in the block /: Ufj = 0,L ^ ^fii^-ii) = Uq(2I + 1). 

(Note that this does not correspond to the total polarization of all the nuclei. The latter 
state would have corresponded to the nuclear density matrix with all blocks I empty except 
for the one with I = N with the only populated state with the grates magnitude of the 
projection L = ±JV). When the pumping is performed with linearly polarized (or non- 
polarized) pulses (i.e., = 7r/4) the steady-state density matrix, as can be shown by direct 
calculations, remains high-temperature Uli — Uq. In the case r ^ oo and complete optical 
orientation (for example when = 0) the Eqs. (37) - (39) allows one to rapidly calculate, in 
a recurrent way, the steady state nuclear density matrix. For that one may use the following 
steps: 

1. using (39) calculate U^ n = 1. 

2. using Eq. (37), successively calculate C/£j, L — 1 — I — 1. 

3. using Eq.(38), derive Uf\ 

Note that the above procedure of calculation of the steady-state nuclear density matrix 
and 300 points of the magnetic-field dependence of the steady-state electron spin polariza- 
tion, for example for a system comprised of 1600 nuclei, for the case of repetitive optical 
orientation of the electron spin takes, when using an up-to-date computer, just a few sec- 
onds. The dynamics can be calculated even faster, because the steady-state density matrix 
is the same for all the time dependence and may be calculated only once. In the numerical 
calculations, one should keep in mind the following. Examination of the function T N (I) (29) 
shows that its width and position of its maximum ~ y/N, i.e., Tn(I) is essentially non-zero 
in the region / ~ \/N. For this reason it is useful, in the numerical calculations to find the 
position of the maximum and, when calculating the sums (27), to restrict oneself to region 
/ < (3 + 4)I max . 

Equations (34) -(36) allow one to calculate the pulse-to-pulse dynamics of the density 
matrix (with all the observables dependent on the nuclear degrees of freedom). For instance, 
the nuclear angular momentum projection (I z ) can be calculated using the formula 

N I 

(^ = E r «WE^ ( 42 ) 

7=0 -/ 
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Note here that the maximum nuclear polarization (for the considered action upon the 
system) is obtained for total polarization inside the blocks with a given J, i.e., Uh = 0,L ^ I 
and Ujj = Uq(2I + 1). Therefore the maximum projection of the nuclear moment, in our 
case, is given by the formula 



N 



N 

{h)ma X = U J2 rV(/)(2J + 1)/ « 2\ - (43) 



7=0 ' 71 



C. Examples of calculating the spin dynamics and the steady-state spin 
polarization for the sigh-constant repetitive spin orientation 

Now, we will present some examples that show that the steady-state polarization and 
the electron spin dynamics substantially depend on the way of preparation of the nuclear 
system and considerably differ on those calculated using the high-temperature nuclear den- 
sity matrix. Figure 2(a) shows magnetic-field dependences of the steady-state electron spin 
polarization under repetitive complete optical orientation for different nuclear relaxation 
times r (curves 1-4). The same figure shows the results of calculations for the case of high- 
temperature nuclear density matrix (curve 0). As the nuclear relaxation time increases, the 
curves differ more and more from the case of the high-temperature nuclear density matrix: 
there arise an asymmetry related to the sign of the optical orientation and a dip in small 
magnetic fields. Figure 2(b) shows a family of magnetic field dependences of the nuclear 
polarization (I z )/N corresponding to figure 2(a). The electron spin dynamics after a long 
repetitive orientation, in the presence of the nuclear relaxation, is shown in figure 3. This 
dynamics, as seen from the figure, is quasi-periodic with the period of the optical orientation 
T. The dynamics in figure 3 corresponds to the case when the train of the pulses, after its 
long action upon the system, is switched off and one can observe the spin dynamics infinitely 
long. If the orienting pulses are not interrupted, then, in the described behaviour the spin 
polarization goes ahead of the orienting pulse, which in some cases, can be indeed observed 
experimentally. 



D. The case of a repetitive sign-alternating spin orientation 

Formulas (34) -(36) are the linear transformation of the vector-column Uj of the nuclear 
density matrix from n-th to n + 1-th pump pulse. In the general case, polarization of 
the orienting pulse may change from pulse to pulse, i.e. = n . Then, the above linear 
transformation can be represented in the matrix form: 

U? +1 = T/(0 n )U™ (44) 

The transformation matrix T/ is given by formulas (34) - (36). When the nuclear relax- 
ation is weak, then in the above case of the complete repetitive sign constant orientation, 
one should expect a strong nuclear polarization (inside the blocks with a given total nuclear 
angular momentum) with decreasing signal of the electron spin dynamics because the elec- 
tron spin, in this case, is conserved. This is the reason why the sign alternating orientation 
of the electron spin is often used in the experiments. In simplest case it means that e.g., for 
the even pulse, = and for the odd pulse <fi = ir/2. In this case the nuclear density matrix 
will always change from pulse to pulse, but, after sufficiently long time, the matrix arising 
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after each even pulse (or after each odd pulse) will remain the same. This matrix evidently 
satisfies the condition 

U, = T j (0)T / (tt/2)U / (45) 

and can be calculated in the same way as it was made above for the case of the sign constant 
orientation. When the nuclear relaxation is absent, the formulas for the non- normalized 
steady-state (in the sense of Eq.(45)) nuclear density matrix have the form: 

1 + $1 T 

Ui - = u -tt^, <«> 

Um = m ± ± ^^iC' = + » L '» L ■ L=l- W -2(47) 

C/, = (7,_i(l + *j) (48) 

The recurrent procedure and normalization are completely similar to those described above. 
For this regime of the orientation, the spin dynamics also exhibits a quasi-periodic behaviour 
(fig. 4a), whereas the magnetic field dependence of the steady-state electron spin polarization 
is getting symmetrical but may show a dip in region of low magnetic fields, (fig. 4b) The 
curve in fig.4b with no dip corresponds to the case of the high-temperature nuclear density 
matrix. 



V. THE ELECTRON-NUCLEAR SPIN ECHO 

As seen from the aforesaid, the system under study is equivalent to an ensemble of 
two-level systems. Each two-dimensional Li-block of the Hamiltonian may be put into 
correspondence with a two-level system with the energy gap Qli- The contribution of each 
LL-block to the electron spin ^-projection is additive and is determined by formula (27). On 
the other hand, this contribution is related to the z-component of the Bloch vector 9% (i) 
(19). For these reasons the system under study should show all the effects inherent of an 
ensemble of two-level systems in the presence of inhomogeneous broadening. 

For example, the effect of arising replicas in spin dynamics (fig. 3), mentioned above, has 
much in common with the effects described in 

QUE Ell- 

Consider now the effect of spin 



echo [13J in such a system. Below, we present a standard description of this effect keeping in 
mind a possible experiment with a quantum dot. Let the electron spin, at t — 0, is optically 
polarized along z-direction in zero magnetic field. Consider the motion of the Bloch vector 
corresponding to some Li-block. Assume, for simplicity, that An = 0. Then the motion 
of the Bloch vector 9lj of the Li-block (11) represents rotation around the rr-axis with the 
frequency flu- 

The Bloch vectors pertaining to different LL-blocks are rotated with different frequencies, 
so that the vectors 9li initially oriented along z-axis, after sufficiently long time span, will 
spread more or less uniformly over the plane perpendicular to the x-axis. [TtJ After it 
happened, let us apply to the system a 7r-pulse of the magnetic field directed along the 
z-axis. In this pulse the magnitude of the magnetic field uj and its duration AT should 
meet the condition wAT = 7r. Denote the initial moment of the 7r-pulse by T\. We assume, 



for simplicity, that the magnetic field is so strong that uj ^> 2Aj_y 1(1 + 1) — L(L — 1) and 
that the dynamics of the Bloch vectors of the Li-blocks, during the pulse, is predominantly 
determined by the magnetic field and, therefore, represents rotation around the z-axis. Thus 
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during interval, AT all the Bloch vectors will rotate around the z-axis by the angle it and 
the vector with the greatest frequency Qli which is, at the beginning of the 7r-pulse, ahead 
of all other vectors, after the pulse will be behind all other vectors. On the contrary, the 
vector with the minimum frequency VLli which is, before the pulse, behind all other vectors, 
after the pulse will be ahead of them. After the end of the 7r-pulse the fastest Bloch vector, 
which is now behind all others, will start to reach the slowest vector which is now ahead 
of others and after the time interval T\ after the pulse they will meet, both being oriented, 
at this moment, along z-axis, as at t = 0. The Bloch vectors between the fastest and the 
slowest ones will also come together after the time interval Ti after the end of the pulse. 
Thus at t = 2Ti + AT, the electron spin projection will be equal to that at t — This is 
what is referred to as the echo signal. The simplifications adopted above do not qualitatively 
change the picture of the effect. It is also clear that instead of the 7r-pulse one can use a 
pulse corresponding to odd number of half-rounds of the electron spin. Formally, the echo 
signal is calculated using formulas (25) and (27), with the needed Bloch vectors of the LI 
blocks calculated using the formulas: 



Here the initial Bloch vector of the Li-block 6 LI (0) is given by Eq.(24) and, thus, depends 
on the nuclear density matrix formed to the moment of observation of the effect. The first 
and the second arguments of the matrices determine there dependence on the magnetic field 
and time, respectively. Figure 5 shows the calculated effect of the 7r-pulse on the electron 
spin dynamics at different values of the pulse amplitude. One can see that the quality of 
the echo signal increases with increasing field amplitude because the dynamics of the system 
approaches the one described above. As is seen from the above treatment, the observation of 
the echo effect, in this case, does not require any inhomogeneous broadening and is possible 
for a single quantum dot. The main difference between the technic described above and the 
standard method of spin echo is that an optical pulse is used instead of the 7r/2-pulse. 
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Captures 

Fig. 1 Topography of the matrix elements of the block with the total angular momentum 
1 = 1. Non-zero matrix elements of the Hamiltonian are colored gray; the blocks with given 
projection of the total nuclear momentum are separated by thick lines. 

Fig.2 (a) Magnetic-field dependences of the steady-state electron spin polarization under 
repetitive complete optical orientation (0 = 0) for different nuclear relaxation times (curves 
1 -4). Curve (0) is calculated for the case of high-temperature nuclear density matrix, (b) 
- the corresponding family of field dependences of the nuclear polarization (I z )/N. 

Fig.3 The electron spin dynamics after a long repetitive orientation of the elctron spin in 
the presence of nuclear relaxation. The arising replicas with the period equal to the orienting 
pulse repetition period is related to the arising odering in the nuclear density matrix. 

Fig.4 Spin dynamics (a) and field dependence of the steady-state spin polarization (b) for 
the case of sign alternating optical orientation and in the absence of the nuclear relaxation. 

Fig. 5 The echo signal in the electron spin dynamics for different values of the 7r-pulse 
amplitude. The step-wise plots show time dependence of the magnetic field. 
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